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Abstract 

The vacuum dynamics of SU(2) lattice gauge theory is studied by means of a 
gauge-invariant effective action defined using the lattice Schrodinger functional. Nu- 
merical simulations are performed both at zero and finite temperature. The vacuum 
is probed using an external constant Abelian chromomagnetic field. The results sug- 
gest that at zero temperature the external field is screened in the continuum limit. 
On the other hand at finite temperature it seems that confinement is restored by 
increasing the strength of the applied field. 
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I. INTRODUCTION 



It is widely recognized that the effective action is a useful tool to investigate the quantum 
properties of field theories. 

In the case of gauge theories when including the quantum fluctuations one faces the 
problem to retain in the effective action the gauge invariance that is manifest at the clas- 
sical level. In the perturbative approach, however, the problem of the gauge invariance of 
the effective action is not so compelling. Indeed in order to perform the perturbative cal- 
culations we need to fix the gauge so that the gauge invariance is lost anyway. Obviously 
the physical quantities turn out to be gauge invariant. In the case of the perturbative 
evaluation of the effective action the problem of the gauge invariance can be efficiently 
resolved by the so-called method of the background effective action In the back- 

ground field approach one separates the quantum field into the fluctuations rj(x) and the 
background field A®(x). In order to define the background field effective action we intro- 
duce the partition function by coupling the external current to the fluctuations. Using 
the background field gauge fixing it is easy to see that the partition function is invari- 
ant against gauge transformation of the background field. In this way, after performing 
the usual Legendre transformation, one obtains an effective action which is invariant for 
background field gauge transformations. 

The lattice approach to gauge theories allows the non perturbative study of gauge 
systems without loosing the gauge invariance. Thus, it is natural to seek for a lattice 
definition of the effective action. Previous attempts (both in three |4],|IJ and four |>] Q 
dimensions) in this direction introduced the background field by means of an external 
current coupled to the lattice gauge field. It turns out, however, that the current term 
added to the lattice gauge action in general is not invariant under the local gauge transfor- 
mations belonging to the gauge group. For instance, if one considers abelian background 
fields then the action with the current term turns out to be invariant only for an abelian 
subgroup of the gauge group. So that only in the case of abelian U(l) gauge theory the 
latttice background field action is gauge invariant. 

The aim of the present paper is to discuss in details a recently proposed method |§ to 
define on the lattice the gauge invariant effective action by using the so-called Schrodinger 



functional gO-pl 



Let us consider the continuum Euclidean Schrodinger functional in Yang- Mills theories 
without matter field: 

Z[AW,A®] = (AW \e~ HT V\ A«) . (1.1) 

In Eq. ( |1 . 1| ) H is the pure gauge Yang-Mills Hamiltonian in the fixed-time temporal gauge, 
T is the Euclidean time extension, while V projects onto the physical states. A^ (x) and 
A^(x) are static classical gauge fields, and the state \A) is such that 

(A|*> = V[A] . (1.2) 

From Eq. Ql.lj) , inserting an orthonormal basis |^ n ) of gauge invariant energy eigenstates, 
it follows: 

Z[A {f \ A {l) ] = J2 e ~ EnT ^n[A {f) ]y*[A {l) ] . (1.3) 
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Note that we are interested in the lattice version of the Schrodinger functional, so that 
it makes sense to perform a discrete sum in Eq. ( |1.3| ) for the spectrum is discrete in a 
finite volume. Eq. ( |1.3| ) shows that the Schrodinger functional is invariant under arbitrary 
gauge transformations of the fields A^> and 

Using standard formal manipulations and the gauge invariance of the Schrodinger 
functional it is easy to rewrite as a functional integral |lO|,p!T|| 



ZMA^JVA^"-™ (1.4) 



with the constraints: 



A^x = 0) = Af 
AJx = T)= 4 f) 



1.5) 



Strictly speaking we should include in Eq. (|1.4j ) the sum over topological inequivalent 
classes. However, it turns out that [12|] on the lattice such an average is not needed because 
the functional integral Eq. flLJ ) is already invariant under arbitrary gauge transformations 
of Af and A ( p. 

On the lattice the natural relation between the continuum gauge fields and the corre- 
sponding lattice links is given by 

Uft, = P exp | iag J dt A^(x + atfi) \ (1-6) 

where P is the path-ordering operator, a is the lattice spacing and g the gauge coupling 
constant. 

The lattice implementation of the Schrodinger functional, Eq. (|1.4j ), is now straight- 
forward: 

Z[U {f \U {i) ] = f VUe~ s . (1.7) 



In Eq. ( |1.7|) the functional integration is done over the links U^x) with the fixed boundary 
values: 

U(x)\ X4=0 = U {{) , U(x)\ X4=T = t/ (/) . (1.8) 
Moreover the action 5* is the standard Wilson action modified to take into account the 



boundaries at x^ = 0,T [12 



^ X,)J,>U 

where U^ u (x) are the plaquettes in the (/i, z/)-plane and 



. 1/2 spatial plaquettes at x 4 = 0, T ,_ N 
W ^ = \ 1 otherwise " ( L1 °) 
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Moreover, it is possible to improve the lattice action S by modifying the weights W^s [|T^ . 
Note that, due to the fact that U® + U^ f \ one cannot impose periodic boundary con- 
ditions in the Euclidean time direction. On the other hand one can assume periodic 
boundary conditions in the spatial directions. 

Let us consider, now, a static external background field A ext (x) = A'^ t (x)X a /2, where 
A a /2 are the generators of the SU(N) Lie algebra. We introduce, now, a new functional: 

nA^] = ~ln{^l\ , (1.11) 



T [ 2(0) 
where 

Z[U cxt ] = Z[U ext ,U cxt ], (1.12) 



and Z[0] means the Schrodinger functional Eq. ( |1.12|) without external background field 



ext. 



(U° xt = 1). The lattice link U ext is obtained from the continuum background field A 
through Eq. flL~6|) . 

From the previous discussion it is clear that rL4 cxt ] is invariant for lattice gauge 
transformations of the external links U^ xt . Morever, from Eq. (|1.3| ) it follows that 



r 



lim T [l cxt ] = E [A cxt ] -Eq[Q] (1.13) 



where E'ofA 6 ^] is the vacuuum energy in presence of the external background field. In 
other words TL4 cxt ] is the lattice gauge- invariant effective action for the static background 
field A ext . In particular, if we consider background fields that give rise to constant field 
strength, then due to the gauge invariance it is easy to show that TL4 ext ] is proportional 
to the spatial volume V. In this case one is interested in the density of the effective action: 



e[l cxt ] = -iln 



3*3 (1 14) 

Z[0] ' 



where Q = V ■ T. We stress that our definition of the lattice effective action uses the 
lattice Schrodinger functional with the same boundary fields at x<± = and = T. So 
that we can glue the two hyperplanes x 4 = and x 4 = T together. This way we end up 
in a lattice with periodic boundary conditions in the time direction too. Therefore our 
lattice Schrodinger functional turns out to be 

Z[U cxt ] = J VU e~ s , (1.15) 

where the functional integral is defined over a four-dimensional hypertorus with the "cold- 
wall" 

U,(x)\ X4=0 = U; xt . (1.16) 



Moreover, due to the lacking of free boundaries, the lattice action in Eq. ( |L15| ) is now the 
familiar Wilson action 

S = S w = \ Tr[l - U^(x)) . (1.17) 
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In this paper we study the properties of the gauge invariant lattice effective action in 
pure gauge non abelian theories. In particular we consider the SU(2) gauge theory in 
presence of constant Abelian chromomagnetic field both at zero and finite temperature. 
The plan of the paper is as follows. In Sect. II we consider the SU(2) gauge theory on 
the lattice in presence of constant Abelian chromomagnetic background field. Section III 
is devoted to the discussion of the Nielsen-Olesen instability on the lattice. In Sect. IV 



we present the numerical results of the Monte Carlo simulations at zero temperature [13 
while Sect. V comprises the finite temperature simulations. Finally our conclusions are 
drawn in Sect. VI. 

II. SU(2) IN A CONSTANT ABELIAN CHROMOMAGNETIC FIELD 

In this paper we are interested in the case of a constant Abelian chromomagnetic field. 
Let us consider the SU(2) gauge theory. In the continuum we have: 

Af(x) = A ext (x)6 a , 3 , At\x) = S k>2Xl H . (2.1) 

The external links corresponding to A ext {x) are easily evaluated from Eq. (|1.6|) : 

Uf* (x , 0) = Ul xt (x , 0) = Ui xt (x, 0) = 1 

(2.2) 

U^\x, 0) = cos (22fE!) + ta 3 sin . 

Our Schrodinger functional Z[A ext ] is defined on a lattice with periodic boundary condi- 
tions, so that we impose that: 

U 2 {xi, x 2 , x 3 , x 4 ) = U 2 {xi + L h x 2 , x 3 , x 4 ) , (2.3) 

where L\ is the lattice extension in the x\ direction (in lattice units). As a consequence 
the magnetic field H turns out to be quantized: 

a 2 qH 2tt 

^- = ^n ext , (2.4) 

with n ext integer. 

According to our previous discussion in evaluating the lattice functional integral 
Eq. ( P-.15 ) we impose that the links belonging to the time slice x 4 = are frozen to 



the configuration Eq. (|2.2| ). Moreover we impose also that the links at the spatial bound- 
aries are fixed according to Eq. ( |2.2j ). In the continuum this last condition amounts to 
the usual requirement that the fluctuations over the background fields vanish at infinity. 
An alternative possibility is given by constraining the links belonging to the the time slice 
X4 = and those at the spatial boundaries to the condition 

U 2 {x) = U^(x) , (2.5) 

while the links U^(x) with ^ ^ 1,2 are unconstrained. The main advantage of the condi- 
tion ( |2.5|) resides in the fact that the time-like plaquettes nearest the frozen hypersurface 
X4 = behave symmetrically in the update procedure. Obviously in the thermodynamic 
limit both conditions should agree as the effective action is concerned. 
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III. THE NIELSEN-OLESEN UNSTABILITY ON THE LATTICE 



As it is well known in the continuum the perturbative evaluation of the effective action 
for the constant Abelian chromomagnetic field faces with the problem of Nielsen-Olesen 
unstable modes |f[4 |. Let us briefly discuss the origin of the unstable modes in the con- 
tinuum. 

In order to evaluate the effective action in the continuum one writes 

A*(x)=Al(x) + rt(x), (3.1) 

where A^(x) = S^^XiH and rj°(x) is the quantum fluctuation over the background field. 
In the background gauge 

[8 a % - g ah %(x)] V l(x) = (3.2) 
we rewrite the pure gauge action in the one-loop approximation as: 

S = S class + d 4 x V ;{x) 0${x) . (3.3) 

The one-loop effective action can be obtained by performing the Gaussian integration 
over the quantum fluctuations and including the Faddeev- Popov determinant. However, 
if we solve the eigenvalue equation 

OfMx) = A^(x) (3.4) 
then we find that there are negative eigenvalues: 

K = vl+vl-gH. (3.5) 

As a matter of fact X u < when gH > p\ + p\. If we perform formally the Gaussian 
functional integration in the one-loop approximation then the effective action picks up 
an imaginary part. The point is that in the functional integration over the unstable 
modes one must include the positive quartic term. It turns out that the unstable modes 
behave like a two-dimensional tachyonic charged scalar field. Thus the dynamics of the 
unstable modes resemble the dynamical Higgs mechanism. As a consequence the response 
of the gauge system to the external field turns out to be strong even in the nominally 
perturbative regime [|i~5f| . 

In order to ascertain if the Nielsen-Olesen one-loop instability survives the lattice 
regularization one should evaluate the Schrodinger functional Eq. ( |1.12j ) in the weak 
coupling region. To this end we write the lattice version of Eq. (ft.lp : 

U ll (x)=exp(tagq ll (x)) U; xt , (3.6) 

where the fluctuations q^(x) = q®(x)a a /2 satisfy the boundary condition 

q^,{x)\ X4=0 = 0. (3.7) 

Inserting Eq. ( |3.6j ) into the plaquette 

U^(x) = U„{x)U v {x + fia)Ul(x + Oa)Ul(x) , (3.8) 
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we rewrite the Wilson action 



X [1>V 



1 - -trU^x) 



(3.9) 



in the quadratic approximation as: 



s w = s ext + s (2) 



where 



5 ext = 4^ 

g 2 



1 — cos 



gHa 2 



(3.10) 



(3.11) 



with Q = Li x L 2 x L 3 x L 4 the lattice volume. Note that the external action S ext in the 
naive continuum limit reduces to the classical action 



tt2 

S d = VT — 



(3.12) 



As concerns the quadratic action S^ 2 \ following the method of Ref. [16] a standard cal- 
culation gives 



a 4 ^{[DM*)-DMx)] 2 U^{x)} 

X,/J.>U 

-a 4 J2 {[DM*), D&WUfiix)} 

X,fl>V 

-2a 2 Tr{[g M (x),g,(x)]^f(x)} 

X,fl>U 

-a 3 Tr {([ D ^(x),q^(x)\ 

-[D^(x), qi/ (x)])U^(x)} 
where is the lattice covariant derivative in the external background U^^x) 



D li f{x) = ~[Uf{x)f{x + jla)U j 



extf , 



X 



m] ■ 



Observing that 



with 



G 



/IV 



1 0^)^(1,2) 

2 



cos(^f) foi/) = (l,2) 



(3.13) 



(3.14) 



(3.15) 



(3.16) 



H, 



(jjl,u)^ (1,2) 

"" 1 sin(^) M = {1,2) 



(3.17) 



6 



and integrating by parts, we rewrite the quadratic action as: 

= 2a 4 G^q^DlD^ - 8 fll/ D*D l7 )q l/ (x)] 

X,/1<1/ 

X,fl>U 

+2m 3 r R[q»(x){Dt + D u )q^x)H^{x)} 

X,fl>U 

-Aia 2 TAq u (x)q^x)H^(x)] , (3.18) 



X,fJ,>U 



where 



D*J(x)= - [f(x) - (UJ?) 1 (x - a£) 

x f(x - ajl)U^(x - ajl)] . (3.19) 

Taking into account that 

= (3-20) 

and using Eq. ( |2.2|) we perform the trace over the color indexes. After a rather long but 
otherwise elementary calculation we get 

= SM(q 3 ) + S ( - 2 \q+,q-), (3.21) 



where 



We have 



qUx)=q 1 Jx)±q 2 Jx). (3.22) 



a: 1 



-ql(x)A* A^Kx)) , (3.23) 



2 

X [IV 



where 



AJ(x) = -[f( x + fia)-f(x) 

(X 



A;f(x) = l[f( x )-f(x-jia)}. 



(3.24) 
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Moreover we have: 



s (2 v, q-) = \ EE G ^ 

X flU 

-qi{x)%T>~q~{xyj 

■ 4 

E (qT(x)v$vfqt(x) - qi{x)V~V~ x q^{xyj H 12 

X 

+1T E (qT(xp2Vfqt(x) - qt{x)V 2 V^q 2 {x)) H 12 

X 

+m 2 E {<h(x)qt(x) - ql( x )qt( x )) #12 , 



(3.25) 



with 



2i 



V±f(x) = A ft f(x)±-5 fl2 sm 



agH 



Xi I f(x + (xa) 



(3.26) 



Obviously we need, now, to fix the gauge. To this end we add a gauge fixing term to 
the action and the associated Faddeev-Popov ghost field action. We use the background 
gauge condition: 

E^(^) = 0. (3.27) 
In the Landau gauge the gauge-fixing term in the one-loop approximation is given by 



(3.28) 



Moreover, in the same approximation we get the following Faddeev-Popov contribution 



SP_ P = -InDet 



(3.29) 



can be extracted from 



The lattice version of the continuum operator Q^ v in Eq. 
Eqs. ( |3.22j ), (|3.25|) , and ( |3.28| ). Unlike the continuum case it is not possible to solve in 
closed form the lattice version of the eigenvalues equations Eq. (|3.4j ). However, if we 
neglect the irrelevant terms and keep only the contributions that survive in the naive 
continuum limit a —>■ 0, then we were able to solve the eigenvalue equations and obtain 
the spectrum. In this approximation we replace G^ u with the identity. So that S^ 2 '(q 3 ) 
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does not depend on the background field and we can discard it. Moreover the sum of 
S^(q + ,q~) and S^) simplifies considerably. We get 

+Y E - ^^Og;^)^^^) + h.c. , (3.30) 

where we restricted the quantum fluctuations to the class of function 

q(x u x 2 , x 3 , x 4 ) = sm(p 4 x i )e i(p2X2+P3X3) f(x 1 ) . (3.31) 

Note that the class of functions Eq. ( |3.31| ) is relevant for the constraint Eq. ( |2.2| ). Similar 
results can be obtained with the constraints Eq. (|2.5|). The periodic boundary conditions 
imply that 

f(x 1 + L 1 ) = f{x 1 ) (3.32) 



2tt 

p^ = —n^ jj, = 2,3,4, (3.33) 
Li,,, 



and n M integer. Within the class of functions Eq. ( |3.31| ) we have 

/n a* a 4 • 2 ( agHxi 
-Ox = -A, Ai + — sin 

cr 

4 

2 A, 2 



H — z y^(l - cosp„a) smfa^iy^i) sm(p 2 a) 

■) 



4 2 / agHx\\ 



2 sin 2 ( - ) (1 - cos(p 2 a)) , (3.34) 



2 = sin( — - — ) + — sin( gHa 2 ) cos(gHaxi) 
er 2 er 

- sin(agHxi) sin(p 2 a) — — — - . (3.35) 



By keeping only the relevant terms, the operators 0\ and 2 further simplify as: 

4 2 / agHxi\ 2 
- 0i = -A x Ai + — sin I — - — I sm{agHxi) sm(p 2 a) 

Ct \ j£i / CL 



2 4 



H — g ^(1 - cos(p 2 a)) , (3.36) 



a 



2z .aV?\ 
— r sin ( 



2? 



sin 



2 ' 
,a 2 gH 



— sm{gHa 2 ) 



(3.37) 
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Introducing the complex scalar fields 

(f) = q-(x) , <f)* = q+(x) 

ijj = qZ(x) , il)* = q+(x) 
where q± = ^(qi ± 92), we get 

2: v=3,4 



(3.38) 



+|^f(x)[-0 1 + (!) 2 ]^)+h.c. 
+TE^( ;r )[- i-WW+h.c. 



(3.39) 



It is easy to verify that the contribution to the one-loop effective action due to the fluctu- 
ating fields q^,fi = 3,4 cancels the one due to the Faddeev-Popov determinant. So that 
we are left with the following quadratic action: 



where 



S(%,0) = o A Y J ^{ x )[-Oi{x)+m 2 mx) 

X 

+a 4 ^0*(x)[-d(a;) - m 2 ]<f){x) 

X 



(3.40) 



Let us introduce the operators 

C = Ai-i 
C* = -A\+i 



e ip 2 a _ -y' 



a 



2 . gHaxi 
+ a Sm( ^ } 



2 . gUax x 
+ a Sm( ^ } - 



(3.41) 

(3.42) 
(3.43) 



By keeping the leading terms in the continuum limit it is not too hard to see that 



1 4 2 

1 = -[C*C + CC*] + ^-(l 



cosp^a) 



(3.44) 



!/=3 

and 

[C,C] = ^ S m(— ). 
So that we find that the eigenvalue equation 

- d/AOn) = \f x ( Xl ) 



(3.45) 



(3.46) 
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admits the solutions: 



An .gHa 2 . 
2 4 

+^ XA 1 - cos ^ fl ) ^ = 0, 1, 2, . . . . (3.47) 

Note that the eigenvalues are degenerate. As a matter of fact the order of the degeneracy 
of the Landau levels turns out to be 

Lia L 2 a 2 . .gHa 2 

We have numerically checked that the approximate spectrum Eq. ( 3.47 ) agrees quite well 
with the exact one as long as Li > 32 and for weak magnetic field. 
From Eqs. ( |3.46|) and fl3.47|) we find the following eigenvalues: 



4 2 



A ™ = ^-^(i-cos^a) 



a 

u=3 



2 ^2 9Ha 2 
2 



+ (2n + 3)— sin 2 ^— — , n = 0, 1, 2, . . . , (3.49) 



4 2 



y=3 



+ (2n - 1)4 sin 2 , n = 0, 1, 2, . . . . (3.50) 

cr 2 



It is now evident that the 0-mode with n = is the Nielsen-Olesen mode with eigenvalues 

2 2 
A u = — (1 - cosp 3 a) + — (1 - cosp 4 a) 
a z a z 

2 . 2 gHa 2 

-^sin (— g— ) I 3 - 51 ) 
which is the discretized version of Eq. Q3.5|) . 

IV. MONTE CARLO SIMULATIONS: T = 

Our numerical simulations have been done on a lattice of size L1L2-L3L4 with periodic 
boundary conditions. In order to project onto the ground state according to Eq.( 1.13|) we 



need L4 3> 1. Moreover in order to be close to the continuum limit Eqs. (|3.11|) and (|2.4| ) 
imply also L\ ^> 1. As a consequence we performed the numerical simulations on lattices 
with Li = L 4 = 32. The transverse size of the lattice L± = L 2 = L 3 has been varied from 
L± = 6 up to L± = 32. We are interested in the density of the effective action Eq. ( |1.14j ). 
We face with the problem of computing a partition function which is the exponential of 
an extensive quantity |T7[. To avoid this problem we consider the derivative of e[A ext ] 
with respect to (3 by taking n cxt (i.e. g H) fixed (see Eq. ( p.4|) ). From Eqs.( |1.14|) , (|1.15|) , 
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and ( p..l7|) it follows: 



,r jexti = de[A^} = _ 1 r 1 gg[^j 

[ J a/3 a [z[u cxt ] dp 



dZ[0) 



Z[0] dp 




(4.1) 

where the subscripts on the average indicate the value of the external links at the bound- 
aries. Obviously e[A ext ] can be obtained by a numerical integration in f3 

i ext nl I j nf I r a ext 



e[A ext ,f3}= d(3'e'[A ext ,f3'], (4.2) 
Jo 

where we have taken into account that lim^o e[A ext , j3] = 0. 

It is evident that the contributions to e'fA 6 ^] due to the frozen time-slice at = 
and to the fixed links at the spatial boundaries must be subtracted. In other words, only 
the dynamical links must be taken into account in evaluating e'[A ext ]. We recall that 
Q = L1L2L3L4 is the total number of lattice sites (i.e. the lattice volume) belonging to 
the lattice A. If we denote with Q ext the lattice sites whose links are fixed according to 
Eq. 



^cxt — L1L2L3 

+ (L 4 - \){L x L 2 Lz - (L x - 2)(L 2 - 2)(L 3 - 2)) , (4.3) 
then the volume occupied by the "internal" lattice sites is given by 

a mt = n- n ext . (4.4) 

Accordingly, we define the derivative of the internal energy density £- nt [A ext ] as 

\ xeA,fi>x 




(4.5) 

4 ext 

where A is the ensemble of the internal lattice sites. 

We use the over-relaxed heat-bath algorithm to update the gauge configurations. Sim- 
ulations have been performed by means of the APEIOO/Quadrics computer. Since we are 
measuring a local quantity such as the plaquette, a low statistics (from 1000 up to 5000 
configurations) is required in order to get a good estimation of e' int . 

In Figure 1 we display the derivative of the energy density normalized to the derivative 
of the external energy density: 

4xt = 1 - c °s(^y-) = 1 - cos(^n ext ) (4.6) 
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versus (3 for L\ = L 4 = 32 and 6 < L± < 32. From Figure 1 we see that in the strong 
coupling region (3 < 1 the external background field is completely shielded. Moreover e' int 



display a peak at (3 ~ 2.2 resembling the behavior of the specific heat |18|,|T9[. This is not 
surprising since our previous studies [|20| in U(l) showed that e[ nt behaves like a specific 
heat. 

In the weak coupling region (3 > 3 Fig. 1 shows that the ratio s[ nt / e' cxt stays constant. 
Actually the constant does depend on L± and n ex t- Indeed in Fig. 1 the dependence on 
L± for fixed external magnetic field is evident. On the other hand in Fig. 2 we keep 
L± = 32 fixed and vary n ex t- We see clearly that the weak coupling plateau decreases by 
increasing the external field. In order to extract Ei nt (f3, n ex t) we can numerically integrate 
the data for e' int (P , n ext ) / £' cxt using the trapezoidal rule 

61*09, ^t) = 4xt f £UP ; neXt) d(3> . (4.7) 
JO £ ext 



In Figure 3 we display e- mt ((3, n ext ) obtained from Eq. (f4~7j). The plateau of the derivative 
of the internal energy density in the weak coupling region results in a linear rising term 
in the energy density. For /? > 1 we get 

e-vatiP, ™ext) ^ {3e' ext a(n ext ) . (4.8) 
Moreover for (3 ^> 1 we get also 

{3e[ nt =(3(1- cos ^-n cxt ) ~ l -H 2 . (4.9) 

So that in the weak coupling region 

£ in t(/0, n ext ) ~ a{n cxt )^H 2 , [3 > 1 . (4.10) 

Figure 1 shows that a(n ex t) — 1 for L± ~ 6 — 8 and n ex t = 1- On the other hand a(n ext ) 
decreases by increasing L± or the external background field. This peculiar behavior can 
be compared with the Abelian case where we found that a(n ext ) — 1 independently on 
L± and n ext 00]. Previous theoretical studies [[HJ suggested that due to the presence of 
the Nielsen-Olesen modes the gauge system reacts strongly to the external perturbation 
even in the nominally perturbative regime. It turns out that the Nielsen-Olesen modes 
behave like a (l+l)-dimensional tachionic charged scalar field. The condensation of these 
modes takes place only in the thermodynamic limit. As a consequence the applied external 
background magnetic field is almost completely screened and there is a dramatic reduction 
of the vacuum magnetic energy. Indeed it turns out that in the infinite volume limit the 
perturbative vacuum and the magnetic condensate vacuum are degenerate for vanishing 
gauge coupling. 

On the lattice the Nielsen-Olesen modes display the one-loop instability when X u given 
by Eq. ( |3.51| ) becomes negative. In the approximation adopted in Sect. Ill we find that 
X u gets negative by increasing L± for fixed external field. Thus we can switch on and 
off the one-loop instability by varying L±. This has been also noticed by the Authors of 
Ref. 0. For instance, by using Eqs. ftOj), ( fU3D , and (|33T|) with L x = L 4 = 32 and 



n ext = 1 we find that X u < for L± > 11. 
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Our numerical results in Fig. 1 show that for n ext = 1 and L± < 10 there is no the 
Nielsen-Olesen instability and the gauge system responds weakly to the external pertur- 
bation in the weak coupling region. On the other hand, by increasing L± (see Fig. 1) or 
n ext (see Fig. 2) we increment the lattice Nielsen-Olesen modes. As a consequence we find 
a clear reduction of the vacuum energy density for both the peak values of e' int (P, n ex t) and 
the coefficient a(n ex t) i n Eq. (|4.10|) decreases towards zero in the thermodynamic limit. 



To see this we need to perform the infinite volume extrapolation. We can extract more 
information from our numerical data by expressing them versus 

(4.11) 



^cff 



where 



is the magnetic length and 



<*- - ^ - i^t„ (412) 



= (4.i3) 



is the lattice effective linear size. Indeed we find that the data for e[ nt (P , n ext ) / e' ext at 
the perturbative tail and at the peak for various lattice sizes and values of n ext can be 
expressed as a function of the scaling variable x (defined in Eq. ( |4.11| )): 



°ext 

For the perturbative tail of e[ nt (/3 , n ext ) / £' ex t we keep the value of the ratio at (3 = 5. On 
the other hand, the peak values have been extracted by fitting the values around the peak 
to (see Fig. 4) 



(4.15) 



4xt 02 (/? -/^peak) 2 + 1 ' 

From Eq. ( |4.15| ) we extract the peak value, oi, and the peak position fl pea k- It has been 



found that the data are compatible with the scaling law Eq. ( 4.14|) with a = 1.5 (see 
Fig. 5). It is remarkable that the same power-law arises if we adopt the alternative 
boundary conditions given by Eq. ( |2.5| ). So that we see that both boundary conditions 
Eq. (|2.3| ) or Eq. (2J3) lead to the same thermodynamic limit. 



If we, further, take into account the shifts A/3 of the peak values, that turns out to 
depend only on n ext , we are led to the universal scaling-law 



x~ a mt ^ ; euj = k((3), (4.16) 



'ext 



where (3 = (3 — A/3. Indeed Figure 6 shows that all our numerical data (for all the values 
of an and L e g) can be approximately arranged on the scaling curve n(j3). Remarkably 
we find that the peak in k{(3) is located at (3 C = 2.2209(68) which agrees with the peak 
position of the specific heat extrapolated to the infinite volume limit f3 c = 2.23(2) fl^| . By 
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using Eq. ([4.16 ) we can determine the infinite volume limit of the vacuum energy density 
e int . We have 



lim e int (/3,n ext ,L eff ) 

L c ff — >oo 



"cxt 



dp' k0') 



x lim 



a H 



^cff 



2P 



dp' k(P') 



x lim 

Lgg—^OO 



a H 



L, 



cff 



0, 



(4.17) 



in the whole range of p. This in turn implies that in the continuum limit (L e g — > oo, P —>■ 
oo) the SU(2) vacuum completely screens the external chromomagnetic Abelian field. In 
other words, the continuum vacuum behaves as an Abelian magnetic condensate medium 
in accordance with the dual superconductivity scenario. 



V. MONTE CARLO SIMULATIONS: T + 

We can extend the study of the SU(2) gauge system in an external chromomagnetic 
Abelian field to the case of finite temperature. As it is well known the relevant dynamical 
quantity is the free energy. On the lattice the physical temperature T p h ys is introduced 
by (in units of Kb — 1): 

7f7— = L t -a, (5.1) 

where L t is the linear extension in the time direction L t = L 4 , while the extension on 
the spatial direction should be infinite. In numerical simulations, however, the spatial 
extension would of course be finite. In order to approximate the thermodynamic limit 
one should respect the relation 

L s > L t . (5.2) 
We perform our numerical simulation on 32 3 x L t lattices by imposing 

T s ~ 4 (5 - 3) 

in order to avoid finite volume effects. 

In the case of constant external chromomagnetic field the relevant quantity is the 
density of the free energy 

F[ ^ ]= __ ln 4_i, v = Ll (5.4) 

The pure gauge system undergoes the deconfinement phase transition by increasing the 
temperature. The order parameter for the deconfinement phase transition is the Polyakov 
loop 

Lt 

p = Tr n u ^ ■ ( 5 - 5 ) 

X4=l 
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As a preliminary step we look at the behavior of the temporal Polyakov loop (P) versus the 
external applied field. We start with the SU(2) gauge system at (3 = 2.5 on 32 3 x 5 lattice 
at zero applied external field (i.e. n ext = 0) that is known to be in the deconfined phase of 
finite temperature SU(2). If the external field strength is increased the expectation value 
of the Polyakov loop is driven towards the value at zero temperature (see Fig. 7). Similar 
behavior has been reported by the authors of Ref. |2l| within a different approach. It is 
worthwhile to stress that our result is consistent with the dual superconductor mechanism 
of confinement. 

On the other hand, if we start with the SU(2) gauge system at zero temperature 
in a constant Abelian chromomagnetic background field of fixed strength (n ext = 1) and 
increase the temperature, then we find that the perturbative tail of the /^-derivative of the 
free energy density F( nt (j3 , n ext ) / e' cxt increases with l/L t and tends towards the "classical" 
value F( nt (p,n ext )/e' cxt ~ 1 (see Fig. 8). 

We may conclude, then, that by increasing the temperature there is no screening effect 
in the free energy density confirming that the zero-temperature screening of the external 
field is related to the confinement. Moreover the information of F( nt (/3, n ext )/e' ext at finite 
temperature can be used to get an estimate of the deconfinement temperature T c . In 
Figure 9 we magnify the peak region for various values of L t . We see clearly that the 
pseudocritical coupling (3*{L t ) depends on L t . To determine the pseudocritical couplings 
we parametrize F( nt (fl, L t ) near the peak as 

F' nt (frL t ) _ ai (L t ) 

e' cxt a 2 {L t )\j3-P*{L t )Y + l K ' } 

We restrict the region near j3*{L t ) until the fits Eq. (|5.6|) give a reduced \ 2 °f order 1. 
Having determined (3*{L t ) we estimate the deconfinement temperature as 

T 1 1 

(5.7) 



Aiatt L t -lf((3*(L t )) 

where 

11 iv 51 '™ , te - 



'<0 = U?i>J «*(-ir">- (5 - 8) 

In Eq.( |5.7| ) we take into account that, due to the frozen time slice, the effective extension 
in the time direction is L e ^ — L t — 1. 

In Figure 10 we display T c /Ai att for different temperatures. Following Ref. [3^j we 
perform a linear extrapolation to the continuum of our data for T c /Ai att . We see that our 
estimate of T c /Ai at t in the continuum is in fair agreement with the one available in the 
literature p3| : 

T 

" = 24.38 ±2.18. (5.9) 



Alatt 

VI. CONCLUSIONS 

We have studied the non-perturbative dynamics of the vacuum of SU(2) lattice gauge the- 
ory by means of the gauge-invariant effective action defined using the lattice Schrodinger 
functional. 
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At zero temperature our numerical results indicate that in the continuum limit L e ^ — > 

oo, (3 — > oo we have 

e[A cxt } = 0, (6.1) 

so that the vacuum screens completely the external chromomagnetic Abelian field. In 
other words, the continuum vacuum behaves as an Abelian magnetic condensate medium 
in accordance with the dual superconductivity scenario. In particular we have 

e[A^] ~ EtBt = -F^ (6.2) 

where /i is the vacuum color magnetic permeability. Thus Eq. ( |6.1[ ) implies that [i — > oo in 
the continuum limit. As a consequence by Lorentz invariance the vacuum color dielectric 
constant tends to zero. This in turns implies that the vacuum does not support an isolated 
color charge, i.e. the color confinement. 

The intimate connection between the screening of the external background field and the 
confinement is corroborated by the finite temperature results. Indeed our numerical data 
show that the zero-temperature screening of the external field is removed by increasing 
the temperature. Moreover, at finite temperature it seems that confinement is restored 
by increasing the strength of the external applied field. 

At finite temperature we find that the /3-derivative of the free energy density behaves 
like a specific heat. From the peak position of the /3-derivative of the free energy density 
we obtained an estimation of the critical temperature T c /Ai att that extrapolates in the 
continuum limit to a value consistent with previous determinations in the literature. 

Let us conclude by stressing that our method can be easily extended to the SU(3) 
gauge theory. Moreover we also feel that the lattice gauge invariant effective action could 
be also employed to study different background fields. 
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Figure 1: The /3-derivative of the internal energy density Eq. ( |4.5|) versus (3 at different 
values of the transverse lattice size L i . 
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Figure 2: The /3-derivative of the internal energy density Eq. (|4.5| ) versus j3 for a transverse 
lattice size Lj_ = 32 at different values of applied external field strength. 
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Figure 3: The energy density Eq. ( |4.7| ) versus (3 for transverse lattice size L± = 32 
at different values of applied external field strength. The solid lines are the linear fits 
Eq. (0). 
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Figure 4: The data displayed in Fig. 2 near the peaks of the /^-derivative of the internal 
energy density Eq. (|4.5|) in correspondence of each value of n ext . The solid lines are the 
fits Eq. (pL15D . 
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Figure 5: The lattice data for e' int (f3 , n cxt ) / e' cxt at different values of the transverse lattice 
size L± and external field strength n ext for (3 = 5 (a) and (3 = j3 pea ,k (b), versus the scaling 
variable defined in Eq. ( [4.111 ). The solid lines are the fits Eq. (|4.14{) with a = 1.5. 
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Figure 6: The scaling curve obtained by re-scaling all lattice data for e[ nt ((3, n ex t) Asff) 
according to Eq. (|4.16|) . 
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Figure 7: The absolute value of the Polyakov loop on a 32 3 x 5 lattice at (3 = 2.5 versus 
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Figure 8: The (3- derivative of the free energy density Eq. Q5.4Q F( nt (/3, n cxt )/e' cxt versus 
at different values of the temporal lattice extension L t . 
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Figure 10: Our lattice data for T c /A latt versus the temperature (circles). The full circle 
is the continuum estrapolation of Ref. |23j (see Eq. (|5.9|) ). The solid line is a linear fit to 
our data. 
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